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Abstract 

It is well known that Forward Time and Centred in Space (FTCS) scheme for scalar 
Hyperbolic Conservation Law (HCL) is unconditionally unstable. The main contribution of 
this work to show that FTCS is conditionally stable for HCL. A new approach is used to give 
bounds on the initial data profile by transforming FTCS into two point convex combination 
scheme. Numerical results are given in support of the claim. 
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1 Introduction 


We consider for the simplicity, the linear transport problem, 

d . . d . . 

—u(x,t) + a—u(x,t) = 0,aG 
at ox 

u o(x) = f (x) 


(1) 

( 2 ) 


where u = u(x, t ) is a scalar field transported by flow of constant velocity a. The exact solution 
of linear transport equation is given by 


u(x, t) = uq(x ± at). (3) 

Seemingly simple, the linear transport equation ([I]) has played a crucial role in the development 
of the numerical methods for general hyperbolic conservation laws. It is the closed form solution 
([3]) of ([l]) which pave the way to analyze any new numerical scheme devised for more complex 
non-linear hyperbolic problems for its convergence and stability. Considering the solution ([3]) , it 
is natural to seek a numerical scheme which is stable in the sense that it does not induce numerical 
oscillations. Therefore, the notion of stability of a numerical schemes evolves around the induced 
spurious numerical oscillations. In order to elaborate it, we consider a uniform grid with the 
spatial width h, time step k and denote the discrete mesh point ( Xj,t n ) by Xj = jh , j £ Z and 
t n = nk, k £ N. 

In the seminal work Courant-Friedrichs and Levy shown that for the convergence a difference 
scheme must contains the physical domain of dependence of partial differential equation |Tj. In 
other words, they gave a necessary condition on ratio of the spatial and time discretization step 
for the convergence of the difference scheme known as CFL number. Since then, the notion of 
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CFL condition has been an indispensable tool for defining the stability of numerical schemes. It 
is known that the linear Von-Neumann stability analysis of a numerical scheme for ([I]) deduces 
a stability condition on the CFL number C = |a|A, A = y Moreover, almost all notion of 
non-linear stability such as upwind range condition, monotone stability unm fl5j . positivity 
preserving mm, total variation stability [4] heavily relies on CFL number to devise any new 
stable scheme e.g., TVD schemes in [3}[7„ 14j, essentially non-oscillatory (ENO) schemes [2, .16], 
weighted essentially non-oscillatory schemes [ T8] . Apart from the notion of stability, in a recent 
paper CFL number is exploited for improved approximation by the flux limiters based scheme in 
[8]. In this work, we are interested in the stability of the following FTCS scheme. It is obtained 
by the discretization of {[lj by replacing the time derivative with a forward difference, and the 
space derivative with a centred difference formula i.e., 

“C = ^ Wh - “?-i) ■ M 

Above three point centred FTCS scheme Q seems to be a correct and natural choice as the 
spatial discretization in FTCS does not violets the physical domain of dependence of ([lj given in 
[T]. It is also interesting to note that the FTCS Q and the centered Lax-Wendroff(LxW) scheme 
m shares the same spatial stencil of grid points. Note that for the CFL number C < 1 the 
three point centred LxW scheme is linearly stable Contrary to the expectation, the solution 
obtained by FTCS scheme Q is diverging and and induced oscillations grow exponentially no 
matter how small the time step is compared to the space step. The classical Von-Neumann 
stability analysis also shows that FTCS (Q is unconditionally unstable. Moreover, FTCS does 
not satisfies any of the notion of non-linear stability mentioned above see j9j. 


2 Non-oscillatory condition on FTCS 


The above, unconditional unstability FTCS scheme Q is surprising as it can be observed that 
for smooth initial data, such as sinusoidal wave, the induced oscillations by FTCS does not grow 
immediately. Moreover, the magnitude and occurrence of induced oscillations can be controlled 
by choosing small CFL number C = A|a|, see Figure [2j On the other hand when applied on 
discontinuous data, FTCS introduces strong oscillations immediately see Figure [3^a). These 
observations have been the motivation for the present study on the dependence of induced oscil¬ 
lations by FTCS scheme on data type and CFL number. In the following result, a classification 
on data type is given in terms of ratio of consecutive gradients which is defined as 


A-u" 

A+uj 

A ~u? 


if a > 0, 
if a < 0, 


( 5 ) 


A ±11,™ = iuj± i =F Uj. We assume that the time and space step ratio satisfies the CFL condition 
aA < 1 and denote the sign of x by sgn(x). 


Theorem 2.1. In the solution data region where 0f € Sftcs = (oo, — 1]U - ^ ^ ^ J —— ,oo 

J [2 — s<?n(a)|a|A 

the FTCS scheme 0 is stable and does not introduce numerical oscillations provided CFL 

number aA < 1. 


Proof. Rewrite FTCS Q in the form 


»” +1 = 


aA 


- - (A + u] + A-u]) . 


( 6 ) 
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Let a > 0: FTCS (§) can be written as, 


u? +1 = v2- 


aX ( A 
T 1 A_?C 


+ 1 A .u 1 ] 


which can be rewritten in upwind stencil form as, 


u] +1 = a + u] + f3 + u]_ v 


where the coefficients or =1- 

2 


aX ( A + id 


A_u™ 


+ 1 and /T = — 


a\ (A + u r - 


A_u™ 


( 7 ) 


( 8 ) 


+ 1 . Note that 


a+f3 = 1, thus to ensure non-oscillatory stable approximation such that Uj_i < — u1 j 

by (|8|), it is sufficient that 


a + = 1 - 


aX ( A + u n . 


2 l A 


i + 1 > 0, and 0+ = 


aA / A + u 71 j 
2 1 A_nl? 


+ 1 > 0 . 


(9) 


Inequalities (10) satisfies if, 


A+< 


A “ 


< 


2 — aA 


aA 


and 


A+< 


A 


> - 1 . 


Which on inversion yield non-oscillatory condition for FTCS scheme ([d]) in case of a > 0, 


A _it' 


n < -1 OR 


A^ 


> 


aA 


A + u!? A + itf 2 — aA 


Let a < 0, then (|6]) can be written in upwind stencil form as, 


■u” +1 = a u'j +/3 Uj +1 . 


. _ L aX ( 

where a = 1 + — 1 + 


A-u™ 

A+v% 


and /3 = 


—aA 


1 + 


A _u r - 


(11) ensures for a non-oscillatory approximation provided 

A _uf ' 


aX 

T 


1 + 


A + u™ 


> — 1 and 


—aA 


1 + 


A _|.u 

A 

A^ 


( 10 ) 

( 11 ) 

Similar to above, 


> 0. 


( 12 ) 


Note that aA < 0, therefore inequalities in (12) satisfy if 


A _< 
1 + A+u 


-2 A 

< —- and 1 + -- J - > 0 


~ aX 


or 


( A _uf 


\A+ii 


< 

n j — 


-(2 + aA) 
aA 


and 


A+u 


A.u" 


n — 


A + u™ 


> -1 


Since > 0, thus on inversion compound inequality satisfies if, 

< -1 OR 


A +Y . , ^ A +< . ~aX 


A_u n ~ 


> 


j 


A_itf ~ (2 + aA) 


(13) 


□ 


The non-oscillatory stable region Sftcs given by inequalities (10) is given for a > 0 in 
Figure ([l]) which clearly shows the effect of CFL number on the stability of FTCS. 
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2.1 Non-oscillatory bound for non-linear scalar problem 

Consider the non-linear scalar hyperbolic conservation laws 


u(x, t) t + g(u(x, t)) x = 0, u(x, 0) = u 0 (x). 


(14) 


The characteristic speed g'{u) associated with (14) at interface x- + i of the cell \xj,Xj+ 1 ] can be 
approximated by 


a 


j+ 


1 = J Ilf ifA + «“#0, 

2 I g\ else. 


(15) 


For (14), FTCS can be written as 


=»” - as (*. - ■ 


(16) 


where gf = gi^uf). Which using (15) can be written as 




°3+\ U j +1 - 


(17) 


In the non-sonic region, following result give bounds for non-oscillatory approximation of non¬ 
linear scalar problem (14) by FTCS. 


Theorem 2.2. Away from the sonic point i.e., where a- i x a-,i > 0, the FTCS scheme (11) 

J 2 J ' 2 _ 

is stable and does not introduce numerical oscillations in the solution region 


on 

'i 



( Q:■7-1-12 

oo, - J 

u 


\ a i -12 


e Sftcs — ‘ 




a i- 12 

oo, —-— 

u 


1, V a J+ 12 J 


Xaj+12 oo) 

if g'(u) 

V 

o 





(18) 

2+Ao 3+ i2 > J 

ifg'(u ) 

A 

o 



provided Amax u \g'(u)\ < 1. 


Proof. Analogous to the proof of Theorem (2.1). 


□ 


3 Numerical Results 


3.1 Linear Case 

Consider linear transport equation ([!]) along with the following initial conditions 

uq{x) = sin(7rx), x G [—1 : 1] 

1 if|®|<l/3 


uq(x) = 


0 


else 


G [-1 : 1 ] 


u 0 (x) = 


ex P^ra 

0 


if \x\ < 1 
else 


x E [-2 : 4] 


(19) 

( 20 ) 

( 21 ) 


In Figure [2j numerical results corresponding to initial condition (19) are given to show the 
effect of CFL number on induced oscillations by FTCS. It can be clearly seen that numerical 
oscillations disappear for small CFL number. This supports the result in Theorem 2.1 as C —> 0 
reduce the region of oscillations of 9 to (-1,0) as shown in Figure [l| 

In order to show data region which cause induced oscillations by FTCS, the following hybrid 
approach using first order non-oscillatory upwind scheme is used to show, 
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• FTCS does not introduce oscillations for data region Sftcs : Use upwind scheme if 0” ^ 
(—1, other wise FTCS is used. Results obtained by this approach (using legend 
’FTCSUP’) in Figure |3](b) and Figure |4](a) show that oscillations does not arise is shown 
as FTCSUP. 

• FTCS introduces oscillations only for data region M.\Sftcs : If 0™ £ (—1, ^r") use FTCS, 
other wise first order upwind scheme is used. Results obtained by this approach in Figure 
|3](c) and Figure |4](b) (using legend ’FTUPCS’) show occurrence of spurious oscillations by 
FTUPCS. 


3.2 Nonlinear scalar 

Consider the non-linear Burgers equations 


,u 


with initial conditions 


u 0 (x) = 


u t T ( ^ )x — 0, 


1 if x < 0.5 

0 x > 0.5, x<G [0,1]. \1 


( 22 ) 


(23) 


Uq(x) = 


(24) 


1 if x = 1 
0, x £ [0, 2] \ 1. 

In Figure [5] and [6] numerical results are given corresponding to initial condition (23) and (24) 
respectively. Similar to the linear case, again results by FTCS and FTUPCS are oscillatory 
whereas results by FTCSUP does not show any numerical oscillations. These results confirm 
that FTCS is conditionally stable and introduce the oscillations only in the solution data region 
where 6 G Sftcs \ ^ 


4 Conclusion 

Non-oscillatory stability bounds on initial data profile are given and numerically tested for un¬ 
conditionally unstable (in Von-Neumann sense) FTCS scheme. An extension of this preliminary 
study for other existing schemes for non-linear hyperbolic problem is being carried out for a 
separate work. 
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a\ 


Figure 1: Plot 9 v/s aX, a = 1: oA —> 0 stabilizes the FTCS as it corresponds to reduced 
oscillatory shaded region. 





Figure 2: 
using N = 


Effect of CFL on induced oscillations in the solution corresponding to (19), at T = 4 
80 grid points (a) C = 0.05 (b) C = 0.25 and (c) C = 0.5. 
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(a) 


(b) 


(c) 


Figure 3: Solution corresponding to (20), with IV = 80 grid point at T = 0.1 and C = 0.1, 
(a) FTCS does not introduce oscillations for top and bottom of left and right discontinuity 
respectively (b) No oscillations by FTCSUP, FTCS is applied only in its stability region and (c) 
Induced oscillations FTUPCS, FTCS is applied only in its unstability region. 




(a) 


(b) 


Figure 4: Solution corresponding to (21), at T = 1.0 and C = 0.6 (a) No oscillations by 
FTCSUP (b) Induced spurious oscillations FTUPCS, FTCS is applied only in its unstability 
region. 




(b) 


(c) 


Figure 5: Solution of Burgers equation corresponding to IC (23) using CFL = 0.9, N = 80, 
after 6 time steps at T = 0.05 
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(a) (b) 


Figure 6: Solution of Burgers equation corresponding to IC 
after 3 time steps 


(23) using CFL = 0.8, N = 40, 
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